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ABSTRACT 

Transiting planets around bright stars have allowed the detailed follow-up and characterization of 
exoplanets, such as the study of exoplanetary atmospheres. The Transit Ephemeris Refinement and 
Monitoring Survey (TERMS) is refining the orbits of the known exoplanets to confirm or rule out 
both transit signatures and the presence of additional companions. Here we present results for the 
companion orbiting HD 114762 in an eccentric 84 day orbit. Radial velocity analysis performed on 
19 years of Lick Observatory data constrain the uncertainty in the predicted time of mid-transit to 
~ 5 hours, which is less than the predicted one-half day transit duration. We find no evidence of 
additional companions in this system. New photometric observations with one of our Automated 
Photoelectric Telescopes (APTs) at Fairborn Observatory taken during a revised transit time for 
companion b, along with 23 years of nightly automated observations, allow us to rule out on-time 
central transits to a limit of ~ 0.001 mag. Early or late central transits are ruled out to a limit 
of ~ 0.002 mag, and transits with half the duration of a central transit are ruled out to a limit of 
- 0.003 mag. 

Subject headings: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (HD 114762) 



1. INTRODUCTION 

The radial velocity (RV) method is still the domi- 
nant source of confirmed exoplanets, although the transit 
technique is rapidly producing planet detections. The 
advantage of transiting planets is that they allow ad- 
ditional characterization of the planet, namely the ra- 
dius and hence the density. For planets transiting bright 
stars, atmospheric studies may be undertaken, a s was 
the case for HP 189733b dKnutson et alJ l2009al) and 
HD 149026b (jKnutson et all l2009bl ). However, most 
of the planets detected using transits orbit relatively 
faint host stars, making them inaccessible to such in- 
vestigations. Transit detections are also biased towards 
short periods due to the geometric transit probability 
(|Kane fc von B raun 20081 ) . The purpose of the Transit 
Ephemeris Refinement and Monitoring Survey (TERMS) 
is to ameliorate these biases through photometric mon- 
itoring of know n RV planets at improved times of pre- 
dicted transits (|Kane et alJl2009T ). 

The first planet candidate to be detected with the 
radial velocity technique orbits the star HD 114762, 
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a late F dwarf. The companion was discovered by 
lLatham et"all (|1989fl who reported the minimum mass 
as 0.011 ± 0.001 Mp). This discovery was confirmed 
by ICochran et all (11991D and a search for transits by 
iRobinson et al.1 ( 19901 ) ruled out t ransit depths greater 
than 0.01 magnitudes. lHalel (|1995| ) used high-resolution 
spectroscopy to measure the projected rotational velocity 
v sin i and concluded that the inclination of the compan- 
ion is likely to be low, possibly rendering the companion 
above the limit for deuterium fusion. There is still uncer- 
tainty as to the true inclination of the companion, par- 
ticularly with the large range of spin-orbit misalignments 
which can occur (|Fabrvcky_fc Winn 200j|. Further long- 
term photometry by I Henry et al.l ( 1997t ) has found the 
host star to be photometrically stable in the V band to 
0.001 magnitudes. 

Here we present our complete RV dataset from Lick 
Observatory that has a time baseline of 19 years. The 
Keplerian analyses of these data is used to produce a re- 
vised orbital solution that further refines the period for 
transit studies and rules out the presence of additional 
companions within the system. We also present 23 years 
of high-precision photometry from the T2 0.25 m and the 
T10 0.8 m Automatic Photoelectric Telescopes (APTs) 
at Fairborn Observatory. The most recent data were ac- 
quired with the T10 APT whose la uncertainty of 0.0012 
mag for a single observations is easily sufficient to find 
or to rule out transits with a predicted depth of ~ 1%. 

2. REVISED ORBITAL PARAMETERS 

The RV data were acquired with the 3.0m Shane tele- 
scope and the Hamilton Echelle Spectrograph at Lick 
Observatory. The data comprise 74 measurements that 
range from 1990 March through 2009 February and are 
shown in Table [T] Of th ese, 46 were extract ed from the 
same spectra as used bv iButler et al.l (|2006l ). However, 
the Lick data reduction pipeline undergoes frequent re- 
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Fig. 1. — Radial velocity data of HD 114762 acquired with the 
3.0m Shane telescope at Lick Observatory. Also shown is the best- 
fit Keplerian model for the planet. 

finement and so each reduction is of superior quality to 
previously published measurements extracted from those 
spectra. As data accumulates, the code adaptively re- 
assesses the relative weights of portions of the spectra 
based upon such factors as contamination by weak tel- 
luric lines or bad pixels. The fourth column in Table [1] 
shows the dewar number which was used with the obser- 
vation, the relevance being the different CCD response 
characteristics which can occur between dewars. This is 
accounted for in the Keplerian orbital fitting described 
below. 

We used Spectroscopy Made Easy (jValenti fc Piskunovl 
119961 ) to fit high- resolution Lick spectra of HD 114862, 
applying the wavelength interva l s, lin e data, and 
methodology of IValenti fc Fischerl (|2005l) . We further 
constrained the surface gravity using Yonsei-Y ale (Y 2 ) 
stellar structure models (iDemarque et alj|2004f) and re- 
vised Hipparcos par allaxes (Ivan Leeuwenll2007l) with the 
iterative method of IValenti et al.l (|2009l ). The resulting 
stellar parameters listed in Table [2] are effective temper- 
ature, surface gravity, iron abundance, projected rota- 
tional velocity, mass, and radius. The stellar radius, 
i?* = 1.24_Rq, is crucial for estimating the depth and 
duration of a planetary transit. These properties are 
consistent with a very metal-poor early- G sub-giant. 

We fit a single-planet Keplerian orbital solution to the 
RV da ta using the techniques described in iHoward et al.l 
(2010) and the partiall y linearized, least-squares fitting 
procedure described in IWright fc Howard! (|2009D . This 
was performed both with and without the inclusion of 
a linear trend in order to determine the significance of 
including that free parameter. The parameter uncertain- 
ties are extracted from the sampling distribution of each 
parameter thro ugh a non-parametric bootstrap analysis 
(|Freed man 1981). Table Ufl lists the fit parameters both 
with and without the trend. 

As indicated by the Xred an d the rms scatter of the 
residuals, the Keplerian fit which includes the trend does 
not improve the fit given the addition of the free pa- 
rameter. Thus this trend is unlikely to be due to any 
additional companions wit hin the system. This is con- 
sistent with the findings of iMugrauer et all (|2006l ) who 
exclude companions with masses greater than 66 Mj at 
orbita l radii of 316-2674 AU. Note that Patience et al. 
( 2002j) detect a s tellar compan i on to HD 114762, further 
characterized by I Bowler et al.l (|2009p. at an angular sep- 
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TABLE 1 
Lick Radial Velocities 



Date Radial Velocity Uncertainty Dewar 



(JD - 2440000) 


(ms 1 ) 


(m s 1 ) 




7964.963900 


-185.78 


22.28 


6 


8017.833900 


315.12 


20.41 


6 


8018.840200 


277.92 


17.40 


6 


8019.830500 


289.41 


27.20 


6 


8375.852800 


-600.95 


14.63 


6 


8437.738000 


262.91 


19.82 


6 


8649.091100 


465.42 


17.25 


8 


8670.973032 


569.49 


21.00 


8 


9068.820700 


452.17 


25.78 


8 


9068.842500 


505.42 


30.66 


8 


9096.814400 


518.85 


27.31 


8 


9096.835900 


532.71 


28.77 


8 


9114.776700 


76.31 


21.17 


8 


9174.705300 


573.39 


14.70 


8 


9350.008000 


476.46 


17.36 


8 


9375.011500 


-296.95 


17.90 


8 


9464.786200 


-636.28 


25.35 


8 


9469.781200 


-554.16 


16.03 


8 


9768.932800 


496.52 


25.17 


39 


9768.955078 


513.54 


22.34 


39 


9801.853516 


-626.62 


24.68 


39 


9801.875000 


-602.62 


28.81 


39 


9801.899414 


-597.07 


24.61 


39 


9802.825195 


-535.91 


23.30 


39 


9802.848633 


-556.53 


20.99 


39 


9802.871094 


-585.50 


23.00 


39 


9803.818359 


-572.41 


21.69 


39 


9803.840820 


-600.02 


22.14 


39 


9803.863281 


-578.70 


24.69 


39 


9858.733398 


409.48 


25.35 


39 


9881.718750 


-445.90 


9.60 


39 


9892.701172 


-304.06 


24.17 


39 


9892.723633 


-317.33 


26.76 


39 


9892.746094 


-336.66 


27.65 


39 


9893.698242 


-218.52 


27.55 


39 


9893.721680 


-271.28 


24.85 


39 


9893.744141 


-256.03 


25.36 


39 


9894.697266 


-237.50 


20.05 


39 


9894.719727 


-170.63 


18.06 


39 


9894.742188 


-203.70 


19.28 


39 


9913.755859 


585.09 


9.80 


39 


9914.675781 


608.46 


9.31 


39 


10072.061523 


355.86 


17.98 


39 


10120.954102 


151.07 


23.88 


39 


10120.976562 


132.44 


25.13 


39 


10121.968750 


96.44 


23.06 


39 


10121.990234 


92.78 


24.77 


39 


10124.940430 


-7.39 


23.47 


39 


10124.962891 


1.09 


26.03 


39 


10126.963867 


-144.75 


26.88 


39 


10126.984375 


-96.89 


39.75 


39 


10128.038086 


-257.56 


24.87 


39 


10128.062500 


-167.79 


26.50 


39 


10128.970703 


-245.66 


22.94 


39 


10128.994141 


-212.46 


20.20 


39 


10144.940430 


-335.55 


11.58 


39 


10172.857422 


716.70 


51.27 


39 


10181.819336 


583.34 


14.87 


39 


10187.833008 


566.28 


23.11 


39 


10187.856445 


520.18 


26.32 


39 


10200.783203 


289.36 


23.24 


39 


10200.805664 


284.61 


22.44 


39 


10504.960938 


623.54 


17.04 


39 


11628.850586 


197.24 


13.89 


18 


12033.864258 


521.20 


10.93 


18 


12833.715820 


-71.10 


9.78 


24 


13068.975586 


-392.18 


14.66 


24 


13132.822266 


373.24 


10.72 


24 


13389.043333 


270.11 


12.56 


24 


13544.723507 


501.86 


10.65 


24 


13545.745370 


494.24 


10.29 


24 


14165.945718 


-625.59 


11.30 


24 


14547.955891 


527.27 


10.38 


24 


14865.043218 


576.89 


13.83 


24 
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TABLE 2 

Stellar Properties 



Parameter 


Value 


V 


7.3 


B-V 


0.525 


Distance (pc) 


38.7 ± 1.1 


T c ft (K) 


5673 ± 44 


log 9 


4.135 ± 0.060 


[Fe/H] 


-0.774 ± 0.030 


usini (kms ) 


1.77 ±0.50 


M* (M ) 


0.83 ±0.01 


R* Ota) 


1.24 ±0.05 



TABLE 3 

Keplerian Fit Parameters 



Parameter 



Trend 



No Trend 



P (days) 

T c a (JD - 2440000) 
T p b (JD - 2440000) 

e 

K (ms" 1 ) 
u> (deg) 

dv/dt (ms -1 yr _1 ) 
Mpsini (Mj) 
a (AU) 

^red 

rms (ms -1 ) 



83.9152 ± 0.0028 
15326.595 ± 0.242 
9889.141 ±0.190 
0.3325 ± 0.0048 
612.72 ±3.35 
201.41 ± 1.04 
-3.83 ± 1.91 
10.99 ±0.09 
0.353 ±0.001 
1.29 
27.53 



83.9151 ± 0.0030 
15326.665 ±0.242 
9889.106 ±0.186 
0.3354 ± 0.0048 
612.48 ±3.52 
201.28 ± 1.01 

N/A 
10.98 ±0.09 
0.353 ± 0.001 
1.32 
27.42 



a Time of mid-transit 
b Time of periastron passa, 



aration of 3.3"(180 AU). This companion could feasibly 
cause the observed trend, but the significance level is too 
low for us to claim such a detection. We therefore adopt 
the solution without the trend included, shown in Fig- 
ure [TJ The fit required four additional free parameters 
due to the offsets induced by the changing of the detec- 
tor dewars. The offsets with respect to the observations 
acquired with dewar 6 are —6.81 ± 5.88, 19.31 ± 6.26, 
17.42 ±8.91, and -5.40 ±7.12 ms" 1 for dewars 8, 39, 18, 
and 24 respectively. The uncertainty in the period for 
the presented orbital solution is a factor of 4 improve - 
ment over the previous solution by iButler et al.l (2006) 
and allows a more robust transit search to be carried out. 

3. TRANSIT EPHEMERIS REFINEMENT 

From the stellar and planetary properties listed in Ta- 
bles[2]and[3j we derive a planetary ra dius of R v = 1.11 Rj 
using the methods described in Bodenhe imer et all 
(2003). This results in a predicted transit duration of 
0.49 days and a predicted transit depth of 0.91%. The 
uncertainty in the stellar mass/radius and subsequent 
uncertainty in the planetary mass/radius have a minor 
effect on the estimated transit duration but in no way 
affects the predicted transit mid-poi nts since these ar e 
derived from the orbital parameters (jKane etaDHQOl). 
Based upon the revised orbital parameters, we computed 
a new transit ephe meris, applied in the follow ing section. 

As described bv iKane fc von Braunl (|2008| ). the prob- 
ability of a planetary transit is a strong function of both 
the eccentricity and the argument of periastron and is 
at a maximum when the periastron passage occurs close 
to the star-observer line of sight (ui = 90°). The mea- 
sured periastron argument for HD 114762b leads to an 
orbital configuration which favors a secondary eclipse 



rather than a primary transit. The transit probability 
based upon the parameters of Table [3] is 1.71%. For 
comparison, an equivalent circlular orbit would have a 
transit probability of 1.73%, but the same eccentricity 
with a periastron argument of 90° would have a transit 
probability of 2.60%. 

4. TRANSIT EXCLUSION AND IMPLICATIONS 

As part of our long-term pr ogram to mea sure luminos- 
ity cycles in solar- type stars ()Henrvl[l999f ). we have ac- 
quired 23 consecutive years of photometric observations 
of HD 114762 with our T2 0.25 m and T10 0.80 m auto- 
matic photometric telescopes (APTs) at Fairborn Obser- 
vatory in southern Arizona. The T2 APT used a photo- 
diode detector to make differential measurements of stars 
through Johnson VRI filters. It collected 1989 nightly 
observations of HD 114762 during the 1989 through 2001 
observing seasons. Details concerning the data acquisi- 
tion and red uctio n with this telesco pe can be found in 
iHenrvl (fl995h and lPercv etldl (I20M ). The T10 APT has 
a two-channel photometer with a dichroic filter that sep- 
arates the Stromgren b and y photometric bands so they 
can be measured simultaneously by two EMI 9124QB 
photomultiplier tubes. T10 collected 1137 nightly obser- 
vations of HD 114762 during the 2002 through 2011 ob- 
serving seasons. We have averaged the Stromgren b and 
y band measurements from T10 to create a new (b + y)/2 
"passband" that gives us improved precision. T10 oper- 
ation i s very simila r to our T8 0.80 m APT, described in 
IHenrvl (fl995l fl999l) . 

The T2 AV and T10 A(b+y) /2 observations are shown 
in the left and right panels, respectively, of Figure [5J 
Both data sets have been normalized so that all seasonal 
means agree with the mean of the first year, shown as 
the horizontal line in the two top panels. This removes 
any long-term brightness variability in the comparison 
stars and in HD 114762 itself; corrections were typically 
a few tenths of a mmag. The top panels of Figure [2] 
give for each APT the standard deviation of a single ob- 
servation from the mean of the whole data set: 0.0122 
mag and 0.0012 mag for T2 and T10, respectively. Both 
values are consistent with the precision of a single ob- 
servation from each telescope, establishing the order of 
magnitude better precision of T10 over T2. Thus, we 
find that HD 114762 is constant to high precision on 
both night-to-night and year-to-year timescales. 

The middle panels of Figure [2] plot the two data sets 
with respect to the orbital phase of the companion, with 
phase 0.0 equal to the predicted time of mid-transit, T c , 
calculated from our revised orbital parameters in Table 
G3 A least-squares sine fit to each data set gives semi- 
amplitudes of 0.0011 ± 0.0004 and 0.00004 ± 0.00005 
mag. This confirms to very high precision the lack 
of stellar activity that might otherwise mimic peri- 
odic radial v e locity yariations in HP 114762 (see, e.g., 
iQuelqz et al.l (120011); IPaulson et al.l (|2004 : iBonfils et all 
(20071): iForveille et al.l (12009^ . 

The bottom panels show the portion of the middle pan- 
els around phase 0.0 plotted with an expanded scale on 
the x-axis. The solid curve in each panel shows the pre- 
dicted time (T c ), depth (~ 1%), and duration (0.49 days) 
of a central transit, computed from the new orbital and 
planetary parameters derived above. The horizontal er- 
ror bar below the transit shows the ± la uncertainty in 
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Fig. 2.— Top panels: The 1989 T2 APT AV observations (left) and the 1137 T10 APT A(b + y)/2 observations (right) of HD 114762 
collected during 23 consecutive years between 1989 and 2011. Middle panels: The T2 and T10 observations plotted modulo the revised RV 
period, P (Tabled with phase 0.0 equal to the predicted time of mid-transit, T c (Table |3). The observations place very tight limits on 
photometric variability on the radial velocity period. Bottom panels: The measurements near phase 0.0 plotted on an expanded scale. The 
solid curve shows the predicted depth (~ 1%) and duration (0.49 days) of a central transit computed with our revised orbital parameters 
(Table El. The horizontal error bar under the transit window gives the ± la uncertainty in T c . Our results show that on-time, central 
transits with a depth of 0.001 mag or more do not occur. 



T c . With the present radial velocity data set in this pa- 
per, the uncertainty in T c is almost exactly equal to the 
transit duration. For the T2 observations in the bottom 
left panel of Figure^ the mean of the 12 in-transit obser- 
vations agrees with the mean of the 1977 out-of-transit 
observations to 0.0001 ± 0.0035 mag, where the uncer- 
tainty is calculated by adding in quadrature the uncer- 
tainty of each of the two means. The uncertainty of the 
in-transit mean level dominates the total uncertainty be- 
cause of the relatively fewer observations in transit than 
out. We have thus measured the difference between the 
in-transit and out-of-transit light levels to a precision of 
0.0035 mag. If we assume the need for a 3c event to 
ensure detection, the T2 observations can only rule out 
transits 0.0105 mag (~ 1%) or deeper. Therefore, the 
T2 observations have just enough sensitivity to detect 
the predicted transits but fail to do so. 

For the much more precise TIO observations in the bot- 
tom right panel of Figure [3J the difference between the 
mean of the 26 in-transit observations and the mean of 
the 1111 out-of-transit observations is 0.0001 ± 0.0003 
mag. Therefore, transits with a 3 a depth of 0.0009 mag 
or more should be detectable. We round this to 0.001 
mag and note that this limit applies only to full-duration, 
on-time events. We must make this distinction because of 
the non-random distribution of observations in the pre- 
dicted transit window. We have only one observation 
very early in the first half of the transit window and 25 
observations that span the second half. Twenty-two of 



these 25 points came from a single night (JD 2455326) 
when we monitored HD 114762 at a higher cadence. The 
worst case for the detection of full-duration transits oc- 
curs if they are 1 standard deviation early, i.e., the win- 
dow slides leftward one-half of its width. In this case, we 
loose all of the monitoring observations from the window 
and pick up only five nightly observations. With only j 
the number of observations in the transit window, the 3 a 
transit depth limit doubles to 0.002 mag. Thus, we can 
confidently rule out early or late full-duration transits to 
a limit of 0.002 mag. If the companion does transit as 
described, this is equivalent to ruling out radii of 0.54i?j, 
or a density of 86 gem -3 . For comparison, the density 
of the planet HAT-P-20b, with a mass of 7.246Mj and a 
radius of 0.867i?j, is 13.8 gem" 3 . 

It is also worth considering transits whose trajectory 
lies along a chord of the star approximately half the stel- 
lar diameter. An on-time transit along such a chord 
would contain half the number of observations, all dur- 
ing the second half of the shorter transit. Thus, the 
detectable limit of an on-time transit with \ the full du- 
ration increases from 0.001 mag to ~ 0014 mag. The 
worst case for detection of these shorter transits occurs 
if they are i a early. In this case, the transit window 
contains only one observation, namely, the first observa- 
tion in the transit window in Figure [5] The standard 
deviation of a single observation is 0.0012 mag, so the 
3 a depth limit for this worst case is ~ 0.0036 mag. 
But this worst case is very unlikely and happens only 
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for T c exactly \ a early. A little earlier or later, and 
the transit window begins picking up additional obser- 
vations, and the depth limit becomes tighter. Therefore, 
we can exclude transits with i full duration to a limit 
of ~ 0.003 mag. Similarly to that described above, this 
rules out transits of companions with densities lower than 
47 gem -3 . 

The T10 observations in the lower-right panel illustrate 
the utility of the TERMS approach to finding transits of 
planets with intermediate orbital periods. Of the 26 in- 
transit observations, only 4 of the nightly observations 
taken over 13 observing seasons fell at random within 
the transit window. The other 22 in-transit observations 
were acquired in a single night, 2010 May 10 UT, pre- 
dicted from our updated orbit based on our new radial 
velocities acquired within the TERMS program. 

5. CONCLUSIONS 

These new RV data refine the orbit of the first plane- 
tary candidate discovered to orbit a solar-type star. Our 
analysis of Lick spectra have greatly improved the esti- 
mated parameters of the host star. Including a linear 
trend into the Keplerian orbital fit does not significantly 
improve the rms scatter of the residuals and we thus 
find no evidence of any additional companions in the sys- 
tem. From the Keplerian fit we produce a refined transit 
ephemeris onto which we have folded 23 years of precise 
photometry acquired with the APT at Fairborn Obser- 
vatory. Our photometric observations rule out on-time 
central transits to a limit of ~ 0.001 mag, early or late 



Bakos, G.A., et al., 2011, ApJ, submitted ( arXiv: 1008.3388 ) 
Bodenheimer, P., Laughlin, G., Lin, D.N.C., 2003, ApJ, 592, 555 
Bonfils, X., et al., 2007, A&A, 474, 293 

Bowler, B.P., Liu, M.C., Cushing, M.C., 2009, ApJ, 706, 1114 
Butler, R.P., et al., 2006, ApJ, 646, 505 

Cochran, W.D., Hatzes, A.P., Hancock, T.J., 1991, ApJ, 380, L35 
Demarque, P., Woo., J., Kim, Y., Yi, S.K., 2004, ApJS, 155, 667 
Fabrycky, D.C., Winn, J.N., 2009, ApJ, 696, 1230 
Freedman, D.A., 1981, The Annals of Statistics, 9, 1218 
Forveille, T., et al., 2009, A&A, 493, 645 
Hale, A., 1995, PASP, 107, 22 

Henry, G.W., 1995, in ASP Conf. Ser. 79, Robotic Telescopes: 
Current Capabilities, Present Developments, and Future 
Prospects for Automated Astronomy, ed. G.W. Henry & J. A. 
Eaton (San Francisco: ASP), 145 

Henry, G.W., Baliunas, S.L., Donahue, R.A., Soon, W.H., Saar, 
S.H., 1997, ApJ, 474, 503 

Henry, G.W., 1999, PASP, 111, 845 

Howard, A.W., et al., 2010, ApJ, 721, 1467 

Kane, S.R., von Braun, K., 2008, ApJ, 689, 492 

Kane, S.R., Mahadevan, S., von Braun, K., Laughlin, G., Ciardi, 
D.R., 2009, PASP, 121, 1386 



central transits to a limit of ~ 0.002 mag, and transits 
of | full duration to a limit of ~ 0.003 mag. Transits 
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orbital inclination. 
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